library(pglm)
FP <- read.csv("Your_Directory/JCR_Jooetal_replication/FP.csv")

m1 <- pglm(pinit~ pinit_lag + gwf_mil*right2 + cinc + majpower + atopally + male + leveledu + 
             militarycareer + combat + milcombat + rebel + irregular, data=FP, family=binomial("logit"), model=c('pooling'))

summary(m1)   

# End of Code